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Abstract 

We study the patterns at which the current flow stabilizes in a ID superconducting wire, for 
various experimentally reasonable boundary conditions, for small fixed current densities and tem¬ 
peratures close to Tc- We pay special attention to the possible existence of a stationary regime. 
If the contacts are superconducting, truly stationary or normal regimes do not exist, but can be 
approached as a limit. In the case of weak superconducting contacts, a rich phase diagram is found, 
with several periodic regimes that involve two phase slip centers. For some of these regimes, the 
density of Cooper pairs does not have mirror symmetry. If the contacts are normal, the stationary 
regime is possible. 

PACS numbers: 74.25.Dw, 74.20.De, 74.25.Sv 
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I. INTRODUCTION: THE QUESTIONS WE INTEND TO ANSWER 


When current is driven along a superconducting wire, it can either flow as normal current, 
as supercurrent, or as a combination of both. The various patterns that may be obtained 
in space and time were reviewed long ago,'i>^ and a new approach^ has been raised in recent 
years. The results obtained for ID wires have been extended to the case of stripes." 

One of the possible patterns is periodic in time, and a salient feature of the periodic 
regime is the appearance of phase slips that occur when the superconducting order parameter 
vanishes at some point. Phase slips occur at definite positions, called phase slip centers 
(PSC). While early studies (e.g. Refs, [ij-yl and references therein) were mainly interested in 
wires of effectively infinite length, so that they extended over many PSC and were insensitive 
to the boundary conditions. Ref. 0 focuses on a parameter region in which there are just a 
few PSC, if any. This is the region that most neatly exhibits the qualitative features of the 
current pattern and, with present microfabrication techniques, it becomes an experimentally 
relevant region. 


Figured] shows the phase diagram for the flow patterns found in Ref. 0 for small currents 
and for temperatures close to Tf.. Below the pertinent curve the sample is in the normal state 
(N) and all the current is normal, whereas above the curve there is a stationary regime (S), 
in which normal current and supercurrent are both present, both are functions of position. 


but none of them depends on time. 
The third possibility found in Ref. y 


.fl 


’he stationary regime was previously found in Ref, 
is the periodic regime, in which the normal current 


and the supercurrent are periodic functions of time; however, following a similar approach 
to Ref. it was found^ that for realistic material parameters the periodic regime is strongly 
disfavored in comparison to the stationary state. 

We must therefore confront an intriguing situation: on the one hand, phase slips are 
described in textbooks and their presence has been confirmed experimentally, whereas the 


stationary possibility is practically ignored in the literature; on the other hand. Refs. 


□ 


3 , 1 ? 


s. One reason for this disparity is the small 


], 7 andO. 


We argue that the other reason 


and la suggest that the stationary regime prevai 
effective length of the wires considered in Refs, 
is that these studies assume that the order parameter vanishes at the boundaries. 

The order parameter does vanish at the boundaries if the “banks” at which the current 
is fed contain ferromagnetic impurities. In the present study we will compare this situation 
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FIG. 1: (Color online) Phase diagram in the current density-temperature plane for the current flow 
regimes, if the order parameter vanishes at the extremes of the wire and u = 5.79. “S” denotes 
stationary regime, “N” denotes exclusively normal current, and “NS” denotes a region where both 
regimes are possible. The blue line is the stability limit when either the current or the temperature 
decreases (i.e., the normal state is stable up to this line), the red line is the stability limit when they 
increase (i.e., the stationary state is stable down to this line), and the purple line is the stability 
limit in both directions. For a wire of length 2L, the values of F in this graph have to be multiplied 
by (L/^(0))“^ and those of j by (L/^(0))“^. 

with other possibilities. The banks will be either superconducting or normal. We may then 
consider the limiting cases of very weak superconductors or of normal metals with very short 
de Gennes length, and examine in what sense the stationary regime is approached. 


II. MODEL 

Following Refs. and 0, we use a minimal model. The wire will be perfectly ID and 
uniform, we will assume electroneutrality, and superconductivity will be described by the 
time-dependent Ginzburg-Landau model (TDGL). We choose a gauge with no vector poten¬ 
tial, write T = 1 — T/Tc and denote by (p the electrochemical potential. The unit of length 
will be denoted by Xq , by to the unit of time, by po the unit of voltage, and by jo the unit 
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( 1 ) 


of the current density. As in Ref. l8|, we take 


.y, 


71 h 4:kBTc AaksTc 

XO = ^ 0 , to = —— , V7o = - , Jo = -. 

SksTc Tie 7re4(0) 

Here ^(0) is the coherence length at T = 0, ks is Boltzmann’s constant, e is the electron 
charge, and a is the normal conductivity. 

With this notation, the TDGL equation and Ohm’s law read 


'ipt = 'ipxx + (r - IV’I^ - V’ 


( 2 ) 


and 

ip^ = ulm - j . (3) 

Here Tp is the order parameter, with normalization imposed by Eq. (E]), the subscripts denote 
partial differentiation with respect to the time t and the arc length x along the wire, and u 
is the ratio between the relaxation times of pj and j.— 

Our model neglects Joule heating and thermal fluctuations. Joule heating can be ne¬ 
glected provided that the current density is sufficiently small and/or the wire is sufficiently 
thin and in good thermal contact with a heat bath. Thermal fluctuations are expected to 
be important near stability limits. 

The wire is assumed to extend along —L < x < L. Equations ([2]) and (j3]) are invariant 
under the transformation x —)■ Lx, t —)■ L'^t, ip —)■ L~^'ip, ip —)■ T L~^r and 

j L~^j, since each of the terms in Eqs. (E]) and ([3]) is multiplied by L~^. Therefore, if also 
the boundary conditions are invariant under this transformation, we can limit our study to a 
single value of L, and the solutions for any other value can be obtained by scaling. Following 
Ref. y, we will set L = 1. 


III. BANKS OF SAME MATERIAL AS THE WIRE 

The most common experimental situation is that the banks and the wire are carved from 
the same layer. The banks are located at the extremes of the wire, x = ±L, and are much 
wider than the wire, so that the current density in them is negligible. The banks can therefore 
be treated as being in equilibrium, and we obtain from Eq. (EJ that |'^(±L)| = T^/^. The 
phase of the order parameter has to obey the Josephson relation and we therefore require 
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the boundary conditions 


^(±L,t) =ri/2 


exp 


—i / (p{±L,t')dt' 


( 4 ) 


The electrochemical potential ip is obtained from Eq. ([3]), in which we will set u = 5.79, as 
appropriate for a superconducting material with paramagnetic impurities in the dirty limit. 

Equations [2H1] were solved numerically; numerical details are provided in the Appendix. 
Figure [2] shows the phase diagram that we encounter in the present case. In contrast to Fig.[T], 
the possible regimes are either purely superconducting (SC) or periodic. It is clear that a 
strictly normal regime is incompatible with the condition (jl]), since the order parameter 
cannot vanish close to x = ±L. Also the stationary regime is implausible: normal currents 
would imply dissipation, so that ip{L) < ip{—L). In this case condition (jl]) implies winding 
that increases with time, and can only be released if the order parameter vanishes at some 
point, so that phase slips are expected. 

It should be mentioned that although in the periodic regime there is a voltage drop along 
the wire, this problem is not equivalent to the case of hxed applied voltage studied in, e.g.. 
In the present case ip{—L) —ip{L) is periodic, rather than constant in time. 


. 


Refs. Id and 
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If there is no normal current and if the order parameter has uniform size, then Eqs. ([2]) 
and ([3]) imply that the maximum current density is jc = 2u{Tjc is called the ‘critical 
current density’, and it is usually stated that the wire is in the normal state for j > jc- In 
Fig. |2]we see that the condition |'0(±L)| = T^/^ favors the superconducting regime for small 
currents, but inhibits it for large currents. 

Since the wire cannot be strictly normal, we investigate in what sense the normal regime 
is approached when j rises above jc- Figure [3] shows the size of the order parameter |'0(x, t)| 
for two current densities above jc- We see that the periodic regime persists, but for larger 
currents phase slips are more frequent, and therefore the order parameter has less time to 
recover and remains small. 


IV. BANKS OF A SUPERCONDUCTING MATERIAL THAT IS WEAKER 
THAN THAT OF THE WIRE 

We consider now the case that the wire and the banks are made of different supercon¬ 
ducting materials, so that |'0(±T,t)| ^ F^/^. We write |'0(±L,t)| = rT^/^ and replace the 
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FIG. 2: (Color online) Phase diagram if the contacts to the wire are made of the same supercon¬ 
ducting material as the wire. “SC” stands for exclusively superconducting current, “PI” stands 
for periodic regime with one phase slip center at the middle of the wire, and “PlSC” denotes a 
region where both regimes are possible. The blue line is the stability limit when either the current 
or the temperature decreases, the red line is the stability limit when they increase, and the purple 
line is the stability limit in both directions. The dashed grey curve corresponds to the critical cur¬ 
rent density jc = 2u(rthat would be obtained in a very long wire in which the boundaries 
become irrelevant. 


boundary condition (jl]) with 


i/>(±L,t) = rF^/^exp 



</9(±L, t')dt' 


(5) 


L -'0 J 

We will especially be interested in the case that r is signihcantly less than 1. 

In this case we hnd a surprisingly rich phase diagram. Figure 0] shows the phase diagram 
obtained for r = 0.05. In addition to the regimes SC and PI, found for r = 1, we hnd 
additional regimes that mediate between them. These regimes are periodic with two PSC at 
symmetric positions with respect to the middle of the wire. Sometimes we found that both 
phase slips occur simultaneously and we denote this regime by P2. We also found cases in 
which the two phase slips are separated in time by half a period, and denote this regime by 
P2'. The stabilization time required to pass between the situations P2 and P2' is very long 
(several times 10^to)- We did not investigate the stability boundaries between P2 and P2' 
for r = 0.05, but this will be done for r = 0.2. 
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FIG. 3: (Color online) Density plot of the absolute value of the order parameter, \^j;{x,t)\, for 
current densities above the “critical current density.” In these plots we took F = 1, so that they 
both correspond to the periodic regime. jV’l = 1 at the boundaries, and vanishes at the phase slip 
points, jc = 2.23. 

The upper stability boundary of PI almost coincides with the upper stability boundary 
of the normal regime in the case '0(±L) = 0, and the lower boundary of P2 almost coincides 
with the lower boundary of the stationary regime. 

It should be born in mind that, for given materials, r is a function of temperature. There¬ 
fore, the phase diagrams for fixed r should not be directly interpreted as phase diagrams in 
the current density-temperature plane. 

As r increases, the region occupied by the purely superconducting regime increases, and 
the topology of the phase diagram can change. Figure [5] shows the phase diagram for r = 0.2. 
We note that in this case there is a region where SC and PI are both possible and that the 
stability boundaries of P2 and P2' have been investigated separately. 

For j < 5.5, the transition between P2 and P2' in Fig. [5] is continuous. In this range 
P2 and P2' are mediated by a regime with two PSC, that we denote by P2", in which the 
time lag changes gradually with decreasing F from zero to half a period. Figure [6] shows 
density plots of \‘ip{x,t)\ for j = 4, for F = 3 (P2' regime) and for F = 3.225 (P2" regime). 
We are not aware of previous reports of situations in a uniform superconducting wire such 
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FIG. 4; (Color online) Phase diagram if the order parameter at the contacts is smaller by a factor 
0.05 than what it would be if the banks and the wire were of the same material. “P2” (or “P2'”) 
stands for periodic regime with two phase slip centers, “SC” and “PI” have the same meaning as 
in Fig. [21 and “P1P2” (or “P1P2'”) denotes a region where one or two PSC are possible. The blue 
line is the stability limit when either the current or the temperature decreases, the red line is the 
stability limit when they increase, and the purple lines are the stability limit in both directions. 
The dashed grey curves are the stability lines for r = 0, presented in Fig. [H 

that the order parameter obeys \'4){—L,t)\ = |'^(L,t)|, but stabilizes in a regime in which 

It is hard to judge numerically whether situations as the case F = 3.225 in Fig. [6] truly 
correspond to a distinct P2" regime, or are rather encountered due to insufficient stabilization 
time. We believe they correspond to a distinct regime, because starting from different initial 
order parameters we obtain the same time lag between consecutive phase slips. 

Since Figs. [Hand [5] show no hysteresis in the transition between SC and P2, we expect a 
continuous transition. At hrst sight this seems impossible, since in P2 the order parameter 
approaches zero at the phase slips, whereas in SC it does not. Indeed, near the phase slips a 
small change in the current density that leads to a transition can result in large differences in 
the order parameters, as illustrated in Fig. [71 However, as shown in Fig. [HI the time between 
consecutive phase slips diverges when the SC regime is approached, so that P2 approaches 
SC almost everywhere in the xt plane. The inset in the hgure shows that the same effect 
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FIG. 5: (Color online) Phase diagram for r = 0.2. The colors of the stability boundaries, as well as 
“P2”, “SC”, “PI” and “P1P2”, have the same meaning as in Fig. U] and, likewise, SCPl and P1P2' 
are regions where two regimes are possible. The stability boundaries that involve two regimes that 
both have two PSC are depicted by dashed lines. The framed region 2 < j < 6 and 2.2 < F < 4 is 
shown in an enlarged scale in the inset. Between the red and the blue dashed lines, P2, P2' and 
PI are possible. 


r=3 r=3.225 
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FIG. 6: (Color online) Density plots of {ipix, t)| for r = 0.2 and j = 4. Darker purplish color means 
smaller \ip\. The case F = 3 is in the P2' regime and the case F = 3.225 is in the P2" regime. 
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FIG. 7: (Color online) Difference between the size of the order parameter in the P2 regime and in 
the SC regime. In both cases we have taken F = 6 and r = 0.05. |?/>| 3.5 denotes \'ijj{x,t)\ for j = 3.5 
(in this case the system is in the P2 regime) and IV’b.s denotes \^{x,t)\ for j = 2.5 (in this case 
the system is in the SC regime). 

occurs when the P2—)-SC transition is due to increase of r. 

An experimentally accessible quantity is the dc component of the voltage drop along the 
wire. Figure |9] shows this quantity for the same parameters as in Fig. [HI which involve a 
hysteresis region. 

For the parameters chosen in Figs. [8] and [9] the position of the PSC in the P2 regime are 
almost independent of j. Also in the P2' regime these positions are almost hxed in most of 
the stability range, but in the case of a continuous transition to PI, very near the transition, 
the PSC migrate to the center of the wire, and thus the PI regime is obtained and the period 
is halved. 

We are now in a position to discuss in what sense the P2 regime approaches the stationary 
regime in the limit r —)■ 0. We hrst note from Fig. [10] that in this limit the PSC approach 
the extremes x = ±L of the wire, and that at these PSC the order parameter V'psc is bound 
to the size iV^pscI < 2rF^/^. Therefore V’(±T) —?• 0 for r —)■ 0, as expected, and there are no 
phase slips in the interior of the wire. 

Next, we recall that a peculiar measurable property of the stationary regime is that 
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FIG. 8: Duration of a period as a function of the current density for the P2 regime (upper curve) 
and for the PI regime (lower curve), for P = 6 and r = 0.05. The curves end where these regimes 
stop to exist. The grey dashed line indicates the transition between P2 and SC. The y-axis is 
logarithmic. Inset; Period in the P2 regime as a function of r, for T = 7 and j = 24; for r > 0.34, 
only supercurrent is present. 



FIG. 9: Time average of the voltage drop along the wire, as a funtion of the current density, for 
F = 6 and r = 0.05. The arrows indicate the direction in which the current is varied. The grey 
dashed line shows the voltage that would be obtained if the wire were in the normal state. 
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FIG. 10: (Color online) The black dots (left y-axis) show the distance beween the phase slip 
centers and the middle of the wire in the P2 regime as a function of r, for F = 7 and j = 24. The 
brown squares (right y-axis) show the maximum size of the order parameter at these phase slip 
centers. The grey line is a guide for the eye and the brown straight line has slope 2. The staircase 
appearance is due to the finite discretization of position. 


the potential felt by Cooper pairs, —ipotodaxg'ip/dt (that can be probed by means of SIS 


junctions, as in Ref. 


12|) is independent of position. Quite generally, we can conclude that 


ii X = xi and x = X2 are two points of the wire such that and r/^(a; 2 ,t) are periodic 

and there is no PSC in the segment xi < x < X 2 , then the time averages of dargip/dt at 
X = Xi and at x = 0:2 have to be equal. Otherwise, | arg'^(xi) — arg'^(x 2 )| would endlessly 
increase, the order parameter would become increasingly wound up, and this would force 
phase slips in the segment xi < x < X 2 - Summing up, we conclude that in the P2 regime 
the time average of daxg'ip/dt is uniform in the segment between the PSC, and in the limit 
r —>■ 0 this property applies to the entire wire. 

In contrast to the stationary regime, in the case of P2 the potential felt by Cooper pairs 
is not independent of position for arbitrary times in the segment between the PSC. Figure 
[m shows the potentials felt by Cooper pairs and by normal electrons between two given 
points, as functions of time, during a period. The normal potential peaks shortly after the 
phase slips. 
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FIG. 11: (Color online) Potential differences between the points —0.7L and 0.7L felt by normal 
electrons and by Cooper pairs, as functions of time, for a situation in the P2 regime. The graph 
extends over one period, and the origin of time was set when a phase slip occurs. The parameters 
used were T = 7, j = 24, and r = 0.3. For these parameters the PSC are at ±0.83L. For a wire 
of length 2L, these parameters correspond to a temperature rc[l — 7(^(0)/L)^], current density 
jo(^(0)/L)^, the potential differences equal the value in the graph times ipo{^{0)/L)'^, and time 
equals the value in the graph multiplied by to(-^/^(0))^. 

V. BANKS OF A NORMAL METAL 

In the spirit of a minimal model description, we adopt the de Gennes boundary condition^ 

= T^(±L)/6 , (6) 

where 6 is a length that represents how far Cooper pairs can survive inside the normal 
metal. The limit 6 —)■ oo corresponds to an insulator and 6 —)■ 0 corresponds to a metal with 
magnetic impurities. Condition (|6]) can be justihed microscopically for static situations, and 
we may expect that it is still qualitatively correct for small current densities. Since condition 
(1^ implies that there is no supercurrent at the contacts, it may also be an appropriate 
description of the case in which quasiparticles are injected and withdrawn from the wire by 
means of NIS junctions. The scaling with length introduced at the end of Sec. [ITl has to be 
supplemented with b —)■ Lb. 

Within a numerical scheme in which position along the wire is replaced by a computational 
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grid and V’ 3 ;(±L) are approximated by finite differences, condition ([ 6 ]) becomes = 

b'ip{:hL =p Ax) /{b + Ax), where Ax is the length of a segment in the grid. 

Unlike the previous sections, and since the points x = are not in equilibrium, the 
rate of change of the phases at these points is not dictated by the potential, so that it is 
possible to have (p{—L) ^ ^{L) without phase slips, and the stationary regime is not ruled 
out. 

The stability boundary of the normal regime can be found as in Ref. y, since the evolution 
equations and the boundary conditions are PT-symmetric,- and all the theorems shown in 
Ref. Q apply here as well. The normal regime is unstable if Eqs. (|2]), ([3]) and ([H]) have a 
solution of the form with Re( 7 ) < T. At the bifurcation from N the 

nonlinear terms in Eqs. ([2]) are ([3]) negligible, and we are left with the spectral problem 

fxx + ixjf = - 7 /, /x(±T) = =f/(±T )/6 , (7) 

where 7 is the eigenvalue with the smallest real part. 

Figure [12] is the phase diagram that we found for b = 0.2L. We see that the boundaries 
move to higher temperatures and the Hopf singularity moves to a lower current density, but 
the topology is the same as in the case 6 = 0 , and the stationary regime still exists. 

Figure [131 shows the potential along the wire when it is in the stationary regime, for 6 = 0 
and for 6 = 0.2, at a representative value of T. In the range j < 40, (p{x) is practically 
proportional to the current density. 


VI. CONCLUSIONS 

We have studied the patterns of current flow along a one dimensional superconducting 
wire, for various boundary conditions, for current densities smaller than 160akBTc^‘^{0)/7reL^ 
and temperatures above Tc[l — 12 ^^( 0 )/L^]. 

A possibility that was predicted long ago,- but to our knowledge has never been detected, 
is the existence of a stationary regime, in which all observable fields are independent of time. 
The stationary regime might find applications in cases in which oscillations would induce 
disturbances (e.g., back action from a SQUID). 

Our results clarify to what extent and in what sense the stationary regime can exist when 
the order parameter is not forced to vanish at the extremes of a superconducting wire. If 
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FIG. 12: (Color online) Phase diagram for normal contacts with de Gennes parameter b = 0.2L. S 
stands for the stationary regime and N for the normal regime. The blue line is the stability limit 
when either the current or the temperature decreases, the red line is the stability limit when they 
increase, and the purple line is the stability limit in both directions. The dashed grey curves are 
the stability lines for 6 = 0, presented in Fig. [TJ 



FIG. 13: Potential along the wire in the stationary regime, divided by the current density, for 6 = 0 
and 6 = 0.2, for P = 10. The thickness of the lines covers the variation of (p{x)jj within the range 
0 < j < 40. The dashed line describes ip{x)/j for the normal regime. 
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the contacts are normal and the de Gennes boundary condition is assumed, the stationary 
regime still exists and the phase diagram is qualitatively unchanged, even if the de Gennes 
length is not negligible in comparison to the length of the wire. 

If the contacts are superconducting, but weaker than the wire, the phase diagram changes 
qualitatively, and several new features arise. The role of the normal regime in Refs. and 
7| is inherited by a periodic regime with one phase slip center, and the order parameter 
convenes nonuniformly to zero as the current density increases. Part of the region in which 


Refs. 


and 0 found the stationary regime, becomes now fully superconducting. The fully 
superconducting regime and the periodic regime with one PSG are mediated by one or 
several periodic regimes with two PSG. In the limit that the contacts are made of very poor 
superconductors, the regime with two PSG becomes similar to the stationary state, since 
the PSG approach the extremes of the wire, and the time average of the potential felt by 
Gooper pairs is uniform between them. 

In the regime with two PSG, the two phase slips may be simultaneous, but there may 
also be a time lag between them. When a time lag is present, the solution of the dynamic 
equations breaks the symmetry of the equations, and |'^(—a:, t)| 7 ^ \%lj{x,t)\. 
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Appendix A: Nnmerical solntion of Eqs. (f4|i 

The wire was discretized into a grid of 2^ vertices and time into steps of 3 x 
We started with an arbitrary function for the order parameter if available, we took 

a function obtained from evolution with similar parameters, and otherwise we usually took 

i) = rPi/T 

Derivatives with respect to x were otained as follows: a fast Fourier transformation was 
performed on "(/(x), and the Fourier transformed function was multiplied by {iky, where k is 
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the reciprocal variable of x and p is the order of the desired derivative. We then transformed 
back to x-space and obtained d^ip/dx^. 

The problem with this method is that it reqnires that V’ be periodic with period 2L, and in 
particnlar that Since this is not necessarily the case, we start by evalnating 

the derivative of = %Ij{x) — L)]x/(2L), which does obey 

and snbseqnently add the derivative of ["^(T) — 'ip{—L)]x/{2L). We are still left with the 
problem that the periodic extension of '0^ is not smooth at ±L and we cannot evalnate its 
derivatives at the bonndaries. However, is not reqnired at the bonndaries, and 'ipxi^L) 
enters the problem after mnltiplication by the length L/2® of a segment, so that V’x(iT) 
can be evalnated as a hnite difference. 

The additional stages of the solntion of Eqs. ([2])-(jl]) pose no nnmerical difficnlty: <p(x) 
can be obtained by trapezoidal integration of Eq. ([3]) and the evolntion of '0 can be followed 
by Euler iteration. Likewise, the phases of can be updated by subtraction of the 

product of the time step times ip{±L). 
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